STOCHASTIC SOLUTIONS OF A CLASS OF HIGHER ORDER 

CAUCHY PROBLEMS IN M'^ 



ERKAN NANE 

Abstract. We study solutions of a class of higher order partial differential equa- 
tions in bounded domains. These partial differential equations appeared first time in 
the papers of AUouba and Zheng Baeumer, Meerschaert and Nane |TU], Meer- 
schaert, Nane and Vellaisamy ^7\, and Nane [32]. We express the solutions by 
subordinating a killed Markov process by a hitting time of a stable subordinator 
of index < /3 < 1 , or by the absolute value of a symmetric a-stable process with 
< a < 2, independent of the Markov process. In some special cases we represent 
the solutions by running composition of k independent Brownian motions, called 
fc- iterated Brownian motion for an integer k > 2. We make use of a connection 
between fractional-time diffusions and higher order partial differential equations es- 
tablished first by Allouba and Zheng [4] and later extended in several directions by 
Baeumer, Meerschaert and Nane [TOl . 



1. Introduction and statement of main results 

In recent years, there have been two hnes of study of the stochastic solutions of 
partial differential equations (PDE's): higher order Cauchy problems [3l |2l HI [101 EH 
132] and time fractional Cauchy problems [9l [131 IMl |37] . We will use the equivalence 
of these two types of Cauchy problems on M'^ and on bounded domains with Dirichlet 
boundary conditions to get classical as well as stochastic solutions of a class of higher 
order Cauchy problems that appeared in (TDl [371 [12] . 

In this paper we suppose that Brownian motion has variance 2t (or that time clock 
is twice the speed of a standard Brownian motion). We express the solutions of these 
Cauchy problems using fc- iterated Brownian motions for an integer k >2. There are 
two ways to define /c-iterated Brownian motions. The first one is just to let 

(1-1) m = B^mi\B,i\...i\B,m---\)\)\) 

where Bj^s are independent real- valued Brownian motions all started at 0. In W^, one 
takes Bi to be an M'^-valued Brownian motion with independent components. In this 
case we denote /c-iterated Brownian motion by /^(t). 

To define the second version of /c-iterated Brownian motions, let X~^{t), X~{t) be 
independent one-dimensional Brownian motions, all started at 0. Two-sided Brownian 
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motion is defined to be 

- I X~{-t), t < 0. 
Then the second version of /c-iterated Brownian motion is defined as 

(1.2) J,(t)=Xi(X2(X3(---(X,(t)) ■■■))) 

where X/s are independent real- valued two-sided Brownian motions all started at 0. 
In M'^, one takes Xi to be an M'^-valued two-sided Brownian motion with independent 
components. In this case we denote fc- iterated Brownian motion by J^it) 

For k = 2, both of these processes I2 and J2 were called iterated Brownian motion 
(IBM) and they have been studied by several researchers; see, for example |3l HI [161 
[Til [281 [391 SOI SH m] and references therein. 

Recently, Aurzada and Lifshits [8] studied the small ball probability for Ik, Jk- 
Arcones [B] studied the large deviation for fc-iterated Brownian motions. 

The classical well-known connection of a PDE and a stochastic process is the Brown- 
ian motion and heat equation connection. Let X(t) G M.'^ be Brownian motion started 
at X. Then the function 

u{t,x) = E,[f{X{t))] 
is the unique solution of the Cauchy problem under mild conditions on / 

d 

— n(t,x) =Au{t,x)] u{0,x) = f{x) 

for t > and x G M*^. There is a similar connection for any Markov process where we 
replace A with the generator of the Markov process; see, for example, [3 [121 [IHl 150] . 

AUouba and Zheng [1] and DeBlassie [2Tj obtained the PDE connection of 2-iterated 
Brownian motion. They showed that for Z{t) = l2{t), or J2{t) 



u{t,x)=E,[f{Z{t))] 



solves the Cauchy problem 



:i.3) ^Mt,x) =^I^ + AMt,x); u{0,x) = f{x) 



for t > and x £ M.'^. The non-Markovian property of IBM is reflected by the 
appearance of the initial function /(x) in the PDE. The methods of Allouba and Zheng 
are more general and they allow one to replace A with the generator of a Markov 
process. Another important characteristic of the AUouba-Zheng work is its setup. Let 
X^ be a continuous Markov process started at x and let B{t) be a Brownian motion 
independent of X. They call X^{\B{t)\) a Brownian time process(BTP). Further 
more: 
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(1) They introduced-for the first time-a large class of stochastic processes includ- 
ing kEBTPs that is obtained by taking at random one of the k copies of 
independent Markov process on each excursion interval of By this 
method some important classes of processes are obtained including iterated 
Brownian motion of Burdzy [16] in the case k = 2, BTPs in the case of = 1, 
Markov snake of Le Gall [301 EOl |32], and a process that is intermediate be- 
tween iterated Brownian motion of Burdzy and Markov snake of Le Gall in 
the case of taking limit — > oo. 

(2) Their formulation and approach handle and link to fourth order PDEs not 
just IBMs, but a much larger class (in (1)) containing many other interesting 
new processes. 

In the absence of the initial function in the PDE fll.3p this problem was studied in 



for t > and x G M'^, where < /3 < 1 and is the generator of some continuous 
Markov process Xoit) started at x = 0. Here d^g(t)/dt^ is the Caputo fractional 
derivative in time, which can be defined as the inverse Laplace transform of s^^(s) — 
s^~^g{0), with g{s) = e^'^^g{t)dt the usual Laplace transform. Zaslavsky [55] used 
flL4l) to model Hamiltonian chaos. 

Baeumer and Meerschaert [9] and Meerschaert and Scheffler [35] show that, in the 
case p{t,x) = T{t)f{x) is a bounded continuous semigroup on a Banach space (with 
corresponding process X(t), E^{t) = mf{u : D{u) > t}, D(t) is a stable subordinator 
with index < /3 < 1), the formula 



yields the unique solution to the fractional Cauchy problem (11. 4p . Here g0{t) is the 
smooth density of the stable subordinator -D(l), such that the Laplace transform 



AUouba and Zheng [1] were also the first to establish a connection between their 
class of BTPs and the time half- derivative through their BTP half-derivative generator 
(see their Theorem 0.5), in addition to connecting their BTPs to 4th order PDEs. 
Essentially, Allouba and Zheng [1] show that BTPs are also stochastic solution of 
(11.41) in the case /3 = 1/2 and L^; is a second order elliptic differential operator of 
divergence form. 

Later, Orsingher and Beghin [13 |16] show that u{t,x) = Kx{f{Ik+i{t))) is the 
solution of 



m 




(1.4) 



u(t, x) 




(1.5) 



u{t, x) 



u{t,x)] 



u{0,x) = f{x) 



dx"^ 
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for t > and x G M. 

We will denote the Laplace, Fourier, and Fourier-Laplace transforms (respectively) 
by: 

POO 

Jo 

u(t,k)= / e~'^'''u{t,x)dx] 

u{s,k)= / e"*'''^ / e~''%{t,x)dtdx. 
Jr'' Jo 

Using Fourier-Laplace transform method, Baeumer, Meerschaert, and Nane [TU] 
showed the equivalence of a class of Higher order Cauchy problems and time fractional 
Cauchy problems: Suppose that X(t) = x + Xoit) where Xoit) is a Levy process in 
M*^ starting at zero. If is the generator of the semigroup p(t,x) = Kx[{f{X(t)))], 
then for any / G D{Lrc), the domain of L^, and for any m = 2,3,4, .. ., both the 
Cauchy problem 

(1-6) -^ = Y.fijj;;;^fJ(^) + L>it,xy, uio,x) = f{x) 

and the fractional Cauchy problem 

Ql/m 

(1.7) -^^u{t,x) = L^u{t,x); m(0,x) = /(x), 

have the same unique solution given by 

POO 

u{t, x) = E,{f{X{E'/"^m) = / piit/sy/"^, x)g,/Us) ds. 

Jo 

Considering (11. 5p . and the equivalence of (ll.6p and (II. 7p . it is natural to expect 
the following Theorem. It establishes the PDE connection of /c— iterated Brownian 
motion which extends PDE connection of 2-iterated Brownian motion (IBM) due to 
AUouba and Zheng and DeBlassie [21]. 

Theorem 1.1. Suppose that X{t) = x + Xo{t) e M.'^ where Xo(t) is a Levy process 
starting at zero. If is the generator of the semigroup T{t)f{x) = Kx[{f{X{t)))], 
then for any f G D{Lx), domain of L^, u(t,x) = E^(/(X(|Jfc(t)|))) is the unique 
solution of the Cauchy problems (ll.6p and (ll.7p with m = 2^ . If Z{t) is a two- 
sided Levy process with independent copies of X for positive and negative times then 
u(t,x) = Kx{f{Z{Jk(t)))) is also the unique solution of the Cauchy problems (11.61) 
and (I1.7P with m = 2^ . 



Remark 1.1. Only trivial extensions to special cases of the formulation in AUouba 
and Zheng [1] lead to the fact that the expressions with I^it) and Jkif) in Theorem 
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\1.1\ are interchangeable in the solution expressions of the given PDEs ( see Theorem 
0.1 and its proof in Allouha and Zheng 

Let -D be a domain in M"'. We define the following spaces of functions. 



C{D) 


= {n 


: D ^R: 


u is continuous}; 


C{D) 


= {u 


: D ^R: 


u is uniformly continuous}; 




= 


: D ^R: 


u is j— times continuously differentiable}; 




= {u 


G C^iD) : 


D'^u is uniformly continuous for all 7 < j} 



Thus, if -u G C^{D), then D'^u continuously extends to D for each multi-index 7 with 

l7l<i- 

We define the spaces of functions C°°{D) = nf^^C^{D) and C^{D) = n°t^C^{D). 

Also, let CJ>(D) be the subspace of C^{D) {C^{D)) that consists of 

functions whose j-th order partial derivatives are uniformly Holder continuous with 
exponent a in D. For simplicity, we will write 

with the understanding that < a < 1 whenever this notation is used, unless other- 
wise stated. 

We use Cc{D),Ci{D),Ci'"{D) to denote those functions in C{D),&{D),&'"{D) 
with compact support. 

A subset D of M'^ is an Z-dimensional manifold with boundary if every point of D 
possesses a neighborhood diffeomorphic to an open set in the space if', which is the 
upper half space in M'. Such a diffeomorphism is called a local parametrization of D. 
The boundary of D, denoted by dD, consists of those points that belong to the image 
of the boundary of if' under some local parametrization. If the diffeomorphism and 
its inverse are C^'°' functions, then we write dD G C-^'". 

Since we are working on a bounded domain, the Fourier transform methods in [36] 
are not useful. Instead we will employ Hilbert space methods used in [37]. Hence, 
given a complete orthonormal basis {ipnix)} on L'^{D), we will call 

u{t,n) = / ipn{x)u{t, x)dx; 
Jd 

I* 1*00 I* 

u{s,n) = I ipnix) / e~^^u{t, x)dtdx = / ^jjn{x)u{s, x)dx. 
Jd Jo Jd 

the ipn, and -i/^n-Laplace transforms, respectively. Since {ipn} is a complete orthonor- 
mal basis for L'^{D), we can invert the ■i/'n-transform 

u{t,x) = y^M(t,n)V;„(x) 

n 

for any t > 0, where the sum converges in the sense (e.g., see [481 Proposition 
10.8.27]). 
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Mittag-Leffler function is defined by Ep{z) = YlT=o r(i+/3fc) ' < /3 < 1; see, for 
example, [33] . 

Let /3 e (0, 1), Doo = (0, oo) X D and define 

f d 



9 , , 



<g{x)t^-\gEL'^iD), t>0}. 



Let Td(X) = inf{t > : X(t) ^ D} be the first exit time of the process X 
from D. We denote by {0„,A„,n > 1} the set of eigenvalues A„ and corresponding 
eigenf unctions 0„ of the Laplacian A with Dirichlet boundary conditions: 

A(j)n = -Xn<Pn m D] (j)\QD = 0. 

We will write u G C^{D) to mean that for each fixed t > 0, u{t, ■) E C^{D), and 
u G C^(-Doo) to mean that u G C'^(-Doo) and is bounded. 

Extending the Fourier-Laplace transform method to bounded domains, Meerschaert, 
Nane and Vellaisamy [37] gave a stochastic as well as an analytic solution to fractional 
Cauchy problem (11. 4p in bounded domains: Let < 7 < 1. Let _D be a bounded 
domain with dD G C^'"' , and T£){t) be the killed semigroup of Brownian motion {Xf} 
in D. Let E^^it) be the process inverse to a stable subordinator of index /3 G (0, 1) 
independent of {X{t)}. Let / G D{/X) n C^{D) n C'^{D) for which the eigenfunction 
expansion (of A/) with respect to the complete orthonormal basis : n G N} 
converges uniformly and absolutely. Then the unique (classical) solution of 

(1.8) ■^u{t,x) = Au{t,x)] xeD,t>0] 

u{t,x) = 0, X E dD, t > 0; 
m(0, x) = f{x), X E D 

is given by 

u{t,x) = ¥.,[f{X{EPmii.rD{X)> E^m 

J. noo POO 

= ^1 Tn{l)f{x)gp{tr^/^)l-^'^~Hl= j Tn{{t/lf)f{x)gp{l)dl 







^1-9) = 5^/»0n(x)i?^(-A„t'^). 



n=l 



Remark 1.2. T/ie analytic representation (11. 9p of the solution is due to Agrawal [T] 
m i/ie case D = (0, M) zs an interval in M. For /3 = 1, t/ie study of the Cauchy 
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problem (11 .Sp boils down to studying heat equation in bounded domains with Dirichlet 
boundary conditions which has solution In this case -E'i(— A„t) = e~^^" . This 

is valid under much less requirements on the initial function and the regularity of 
the boundary, see for example |12]. The solution to (II. 8p is also given by (II. 9p if 
we replace Brownian motion with a diffusion process in which case the Laplacian A 
should be replaced with the diffusion operator, see Meerschaert, et. al. |37j . 

Let 

H^ADoo) = !^u:D^^R: \^u{t,x)\ < g{x)t'/^~\ g e L^{D), t > 0, 

^u,A^u G C(D^),k = I,-- - ,m,A'=M G C\D),k = I,-- - ,m - ll. 
ot J 

Using the equivalence of fractional Cauchy problem (II. 8p with the Higher order 
Cauchy problems (II. 6p with the correct Dirichlet type Boundary conditions we obtain 
the second main result in this paper. 

Theorem 1.2. Let m = 2,3, ■ ■ ■ , be an integer. Let D be a bounded domain with 
dD G C-^'^ , and T^lt) be the killed semigroup of Brownian motion {Xf} in D. Let 
{^^/'"(t)} be the process inverse to a stable subordinator of index 1/m independent 
of{X{t)}. Let f G D(A) nC2™-3p) p^2m-2p)(^^ L'^{D)) be such that the eigen- 
function expansion of A"^^^f with respect to : n > 1} converges absolutely and 
uniformly. Then the (classical) solution of 



:i.lO) = ^^^AV(x) + A-«(t,x), xGAt>0; 

dt ^ r(j/m) 



u{t, x) = A'u(t, x) = 0, t > 0, X edD, / = 1, ■ ■ • m - 1; 

m(0, x) = f{x), X E D 

is given by 

u{t,x) = E.[/(X(i?V-(t)))J(r^(X)>Ei/-(t))] 

(1.11) = E4f{X{E'/^it)))HrD{X{E'/"^))>t)] 

POD 

= tm / TD{l)f{x)gi/^{tl-"')l-"'--'dl 
Jo 

CO 

= 5^/(n)0„(x)£;i/^(-A„t^/"^). 

n=l 

In the case m = 2^ for some integer k > 1, the solution to (ll.lOp is also given by 
u{t,x) = E^[f{X{\Ik{t)\))I{TD{X) > \Ik{t)\)]. For (3 = 1/2^= the solution to (OD is 
also given byu{t,x) = E,[/(X(|4(t)|))/(rz5(X) > |4(t)|)]. 
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Meerschaert, Nane and Vellaisamy [37] proved this theorem in the case m = 2. 



Remark 1.3. Theorem \1.2\ also holds with the version of k -iterated Brownian motion 
Jk(t). Here, the outer process X{t) is a two-sided Brownian motion and Jk{t) is an 
independent k-iterated Brownian motion. In this case, using a simple conditioning 
argument, we can show that the function 

u{t,x) = EMXX{Jkmii-rD{X~) < ,h{t) < r^(X+))] 

reduces to Equation fll.lll) and hence is also a solution to both Cauchy problems (11. 101) 
and (11.81) with m = 2^ . 



In [l2], we studied the Cauchy problems that can be solved by running a-time 
processes with < a < 2. An a-time process is a Markov process in which the time 
parameter is replaced with the absolute value of an independent symmetric a-stable 
process Y with < a < 2. 

As a special case in Nane Theorem 2.1] we established: Let {X{t)} be a 
continuous Markov process with generator A, and let {Y{t)} be a Cauchy process 
independent of {X{t)}. Let / be a bounded measurable function in the domain of 
A, with Dijf bounded and Holder continuous for all I < i, j < d. Then u{t,x) = 
E4f{X{\Y{t)\)] is a solution of 

^u{t,x) = -MM _ AMt^x); u{0,x) = f{x) 
at-^ nt 

for t > 0, X G W^. 

This reduces to nonhomogeneous wave equation in the case X is another 
Cauchy process independent of y, the generator A = — (— A)^/^, fractional Laplacian, 
i.e., u{t,x) = Kx[f {X {\Y {t)\)] is a solution of 

^u{t,x) = M^)^!M + Au{t,x); u{0,x) = f{x) 

Ot^ Tit 

for t > 0, X G M'^. This is one of the most interesting PDE connections of these 
iterated processes. 
Let 

/CA2 = {n:D^^M: — n, An, A^n G ^(Z^oo)} 
In bounded domains we obtain the following 

Theorem 1.3. Let D he a hounded domain with dD G C^''^, < 7 < 1. Let {X(t)} he 
a Brownian motion in with independent components, and let {Y{t)} he a Cauchy 
process independent of {X{t)}. Let f G -D(A) fl C'^{D) for which the eigenfunction 
expansion of Af with respect to the complete orthonormal hasis {0„ : n > 1} con- 
verges absolutely and uniformly. Then u(t,x) = E^[/(X(|F(t)|)/(r£)(X) > is 
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a solution of 



(1.12) ^u(t,x) = _ A\(t,x), t>0, xeD 

u(0, x) = f{x), X E D, 

u{t,x) = AM(t,x) = 0, xedD,t>0. 

Remark 1.4. In Theorem the solution is expressed by subordinating a killed 
process by an increasing process E^^'^it), inverse process to a stable subordinator 
of index 1/m. And we see from Theorem \l.^ that killing the outer process and then 
subordinating, or subordinating and then killing gives the same solution of the Cauchy 
problem fll.lOp . The solution in the case m = 2'^ is also given by subordinating the 
killed process by k-iterated Brownian motion. Is there an increasing process A{t) such 
that we get a similar relation for the solution of f ll.l2p ? 

Remark 1.5. We discuss how to extend Theorems and \l.!^ to other contin- 
uous Markov processes in section 6. 

This paper is organized as follows. In section 2 we give some preliminaries. Section 
3 is devoted to the proof of Theorem ll.il We prove Theorem ll.2l in section 4. Theorem 
11.31 is proved in section 5. We also state and prove a theorem for Cauchy problems 
that can be solved by running a Brownian motion subordinated to the absolute value 
of a symmetric a-stable process with a G (0, 2) a rational, and a ^ 1. We discuss 
extensions of the theorems proved in this paper to other types of Markov process in 
section 6. 



2. PRELIMINARIES 

Let Xo(t) be a Levy process started at zero and X{t) = x + X^it) for x G M'^, the 
generator of the semigroup T(t)f{x) = E^.[/(X(t))] is a pseudo-differential opera- 
tor [5l |271 152] that can be explicitly computed by inverting the Levy representation. 
The Levy process Xo{t) has characteristic function 

E[exp{ik ■ Xoit))] = exp{tij{k)) 

with 

^^k) = tk.a-h-Qk+ [ (e^'-y - 1 - j^^irTi?! ^(^2/)' 

where a G M"^, Q is a nonnegative definite matrix, and is a a-finite Borel measure 
on M'^ such that 

/ min{l, ||?/||^}z/((i?/) < oo; 

Jy^O 



see for example [5^ Theorem 3.1.11] and |5l Theorem 1.2.14]. Let 

/>)= / e-''-'f{x)dx 



denote the Fourier transform. Theorem 3.1 in [9] shows that Lxf{x) is the inverse 
Fourier transform of ■ip{k)f{k) for all / G D{Lx), where 

D{L,) = {fe L\R'') : ^{k)f{k) = h{k) 3 he L^R")}, 



+ / (fix + y)- fix) - ) uidy) 



and 



(2.1) 

-4_ / I {(t -I- 7/^ - f('r.'\ - ZJJ^ 

i + y' 

for all / G iy^'^(M'^), the Sobolev space of L^-functions whose first and second partial 
derivatives are all L^-functions. This includes the special case where Xo{t) is an oper- 
ator Levy motion. We can also write = ■?/;(— iV) where V = [d/dxi, . . . , d/dxd)' ■ 
For example, if Xo(t) is spherically symmetric stable then il){k) = — -D||A;||" and Lr^ = 
— A)"/^, a fractional derivative in space, using the correspondence kj — > —id/dxj 
for 1 < j < d. 

We say that D satisfies an exterior cone condition at a fixed point xq G dD if there 
exists a finite right circular cone V = Vxq with vertex xq such that D fl \4o = ^o, and 
a uniform exterior cone condition if D satisfies an exterior cone condition at every 
point Xq G dD and the cones V^p are all congruent to some fixed cone V. 

Let {Xt} be a Brownian motion in and A C M"'. Let Ta = inf{t > : Xt e A} 
be the first hitting time of the set A. We say that a point y is regular for a set A if 
Py[TA = 0] = 1. Note that a point y is regular for a set A provided, starting at y, the 
process does not go a positive length of time before hitting A. 

The right condition for the existence of the solution to the Dirichlet problem turns 
out to be that every point of dD is regular for D*-^ (cf. [121 Section II. 1]). 

If a domain satisfies a uniform exterior cone condition, then every point of dD is 
regular for D'-' . 

Let dD G C^. Then at each point x G dD there exists a unique outward pointing 
unit vector 

e{x) = i9^{x),--- ,ed{x)). 

Let u G C^{D), the set of functions which have continuous extension of the first 
derivative up to the boundary. Let 

du 

denote the directional derivative, where Vm is the gradient vector of u. 

Now we recall Green's first and second identities (see, for example, [Ml Section 
2.4]). Let u,ve C'^{D)nC\D). Then 
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r Ou f f dv 

/ -- — vdx = / uvOids — I u— — dx (integration by parts formula), 
J D dxi Jqd Jd oxi 

[ \/v-\/udx = — [ uAvdx + [ ^uds (Green's first identity), 
Jd Jd JdD 

I [uAv — vAu\dx 
Jd 



dv du 



ds, (Green's second identity). 



IdD 

Let D be bounded and every point of dD be regular for . Markov process 
corresponding to the Dirichlet problem is a killed Brownian motion. We denote the 
eigenvalues and the eigenfunctions of A by {Xn, (pn}'^=i, where e C°°{D). The 
corresponding heat kernel is given by 

CXI 

PD{t,x,y) = ^e"^"Vn(x)(/)„,(y). 

n=l 

The series converges absolutely and uniformly on [to, oo) x x for all to > 0. In 

this case, the semigroup given by 

(2.2) 

TD{t)f{x)=E4f{Xt)I{t<TD{X))]= / pz?(t,a;,y)/(y)rfy = Ve"^"Vn(x)/» 

solves the Heat equation in D with Dirichlet boundary conditions: 

^'^^l'^^ = Au(t,x), xeD,t>0, 
dt ^ ^ 

u{t, x) = 0, X E dD, 
m(0,x) = /(x), X E D. 

Remark 2.1. The eigenfunctions belong to {D) r\ {D) , by pH Corollary 8.11, 
Theorems 8.15 and 8.24]. If D satisfies a uniform exterior cone condition all the 
eigenfunctions belong to C°'{D) by [2U Theorem 8.29]. If dD G C^'", then all the 
eigenfunctions belong to C^'°'{D) by [211 Corollary 8.36]. If dD G then each 
eigenfunction of A is in C°°{D) by Theorem 8.13]. 



Inverse stable subordinators arise in [23 ESI ESj as scaling limits of continuous time 
random walks. Let D{t) be a stable subordinator of index < /3 < 1. The hitting 
time of D is defined as E^(t) = inf{a; : D{x)>t}, so that {^^(t) < x} = {D{x) > t]. 
We will call the process inverse to stable subordinator D in this paper. Writing 
gi3{u) for the density of -D(l), it follows that D{t) has density t~^/^ gj3{t~^/^u) for any 
t > 0. Using the inverse relation P{E(t) < x) = P{D{x) > t) and taking derivatives, 
it follows that E{t) has density 

(2.3) /,(x)=t/3-V-'/%(te-^/^), 
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whose t (— 7- s Laplace transform can also be derived from the equation 



oo 



by taking Laplace transforms on both sides. 

3. Cauchy problems 
We prove Theorem 11.11 in this section. First we need the following Lemmas. 

Lemma 3.1. Let E'^'^,E^^ be two independent processes that are inverses to stable 
subordinators of index < Pi,P2 < 1- Then E^^^ [E^'^ {t)) is inverse to a stable subor- 
dinator of index /3i/32- 

Proof. Since composition of stable subordinators gives another stable subordinator, 
see Bochner [T5], the result follows as E^^ is the inverse of a stable subordinator. □ 



Lemma 3.2. For fixed t > 0, k-iterated Brownian motion 

\h{t)\ = \B,{\B,{\...{\B,m...\)\)\ 
and E^/"^^ [t) have the same one- dimensional distributions. 

Proof. Since ^1/2 and \B{t)\ have the same 1-dimensional distributions, see, for ex- 
ample, proof of Theorem 3.1 in [TD], the proof follows from Lemma [3.11 by induction 
on k. Hence we have by composing k independent E^/'^s 

EnEn..{Ei'\t))))=E^'^\t)'^\h{t)\. 

□ 

Corollary 3.1. 

4+1 (t) = B,{\B,{\B,{\ . . . {B,^,{t)) ■ ■ ■ DDI) 

Proof of Theorem \l.l[ The proof is an adaptation of the proof of Theorem 3.1 in 
Baeumer et. al. [10]. Using Lemma [3^2] we get that 

uit^x) = E.(/(X(|4(t)D)) = E{f{X{Ey'"))) 

is a solution to both the Higher order Cauchy problem fll.7l) and the fractional Cauchy 
problem f 1 1.6 1) for m = 2^. By a simple conditioning argument we also get that 
u{t,x) = E^{f{Z{Jk{t)))). □ 
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4. Cauchy problems in Bounded domains 



The inverse stable subordinators with /3 = 1/2^^ are related to Brownian subordi- 
nators by Lemma [3.2[ this is well-known for the case k = 1, see, for example, [10]. 



Since Brownian subordinators are related to higher-order Cauchy problems by The- 
orem II. H this relationship can also be used to connect those higher-order Cauchy 
problems in bounded domains to their time-fractional analogues. In this section, we 
establish those connections for Cauchy problems on bounded domains in M.'^. We 
extend this to establish an equivalence between a killed Markov process subordinated 
to an inverse stable subordinator with f3 = 1/2^', and the same process subject to 
a Brownian subordinator in section 6. Finally, we identify the boundary conditions 
that make the two formulations identical. This solves an open problem in [10]. This 
problem was solved in [STJ Theorem 4.1] for k = 1. 



Lemma 4.1. Let D be a bounded domain with dD G C^'", < a < 1. Let {</>„, A„, n > 
1} be the set of eigenvalues A„, and corresponding eigenf unctions of the Laplacian 
A. Let f G C^{D) for j = 1,- ■ ■2k and all the partial derivatives of f of order up to 
2k — 1 vanish on the boundary (A simpler condition is to assume f G Cl^{D) ). Then 

l,---fc 



(4.1) / 0„(x)AV(x)rfx = (-A„)^ / 0„(x)/(x)rfx= (-A„)V(ri), J 
J D Jd 

Proof. We use Green's second identity and induction in j 

dA^-^f 



6„AV - A^-'fA(l)^]d. 



■X 



de 



89 



ds. 



where we use the fact that A^ f\gE, = = 4>n\dD, f £ C^{D) for j = 1, ■ ■ ■ ,2k, 
and (pn G C^''^{D) by Remark 12.11 Hence, by induction, the 0„-transform of A^u is 
{—^nY f{n), as (pn is the eigenfunction of the Laplacian corresponding to eigenvalue 
A„. □ 



Proof of Theorem \l.S[ Suppose m is a solution to Equation ( ll.lOp . Taking the 
transform of f ll.lOp and using Lemma 14. H we obtain 



(4.2) 



d 



m—l 



p/m-1 



—u(t,n) = > ^, , . 
dt > ^ r(j/m) 



-KyKn) + {-\nru{t,n). 



Note that the time derivative commutes with the 0„-transform, as 



Q-n{t,x) 



< g{x)t^^'"~\ g e L'^iD), t > 0. 



Taking Laplace transforms on both sides and using the well-known Laplace trans- 
form formula 

r-oo ^-fS 



(4.3) 



r(i-/3) 



e-'^dt = s 



0-1 
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for /3 < 1, gives us 

m— 1 

(4.4) su{s,n) - u{0,n) = ^ s-^/'"(-A„)V» + (-A„)'"m(s, n). 

i=i 

Since u is uniformly continuous on C([0,e] x D), it is also uniformly bounded on 
[0, e] X D. So, we have limt^o Jiju{t,x)(f)n{x)dx = f{n). Hence, u(0,n) = f{n). By 
collecting the like terms, we obtain 

Rn) (l + Er=i'^"''/™(-An)0 

(4.5) u{s,n) = ^ —— ^. 

s- (-A„)- 

Using the simple equality (a - 6)(a™-^ + + ■ + ab^^^ + b^"^) = a™ - 6™ for 

any m = 2, 3, ■ ■ ■ , for fixed n and for large s, we get 



(4.6) u{s,n) 



s- (-A„)™ 



(4.7) 



(sV™ - (-A„)) (si-V- + ^--1 
si/'"-7(n) 



sV™ + A„ ■ 

It follows from the proof of Theorem 3.1 and Corollary 3.2 in [37] that the inverse 
0„-Laplace transform of (14.71) is given by 

oo 

(4.8) u{t,x) = 5^/»^i/™(-A„ti/™)</.„(a:) 

n.=l 

= E4fiX{E'/"^mi{rn{X)>E'/"^m 
= E.[/(X(i?i/-(t)))/(r^(X(i?i/-)) >t)], 

where E^f'^(t) is the process inverse to stable subordinator of index 1/m. 

For any fixed n > 1, the two formulae (14.51) and (14. 7p are well-defined and equal 
for all sufficiently large s. Since the inverse Laplace transform of u{s,n) in (14. 7p is 

(4.9) u{t,n) = fin)E,/U-Kt'^n, 

see, for example, [33], we can see easily that u{t,n) is continuous in t > for any 
n > 1. Hence, the uniqueness theorem for Laplace transforms [TJ Theorem L7.3] 
shows that, for each n > 1, u{t,n) is the unique continuous function whose Laplace 
transform is given by (14.61) . Since x i— )■ u(t,x) is an element of L'^{D) for every 
t > 0, and two elements of L^{D) with the same 0„-transform are equal dx-almost 
everywhere, we have (14. 8 p is the unique element of L'^{D) and (14. 9 p is its 0„-transform. 
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m. 



Now, we show that the solution u defined by f l4.8|] satisfies all the properties in 
fll.lOp . From the Proof of Theorem 3.1 in [37], we can get that the solution u{t,x) 
defined by the series (14. 8 p belongs to L^{D), converges absolutely and uniformly for 
^ > ^0 > for some to- 
Next we show that A'm G C(-Doo) for / = 1, ■ ■ ■ ,m. To do this, we need only to 
show the absolute and uniform convergence of the series defining A'm for / = 1, 
To apply A' term-by-term to (14. Sp . we have to show that the series 

oo 
n=l 

is absolutely and uniformly convergent for t > to > 0. 

Note that by Lemma [4. II the 0„-transform of f is given by 

0„(x)A-'/(x)rfx = {-XnYfin) 

D 

and using [291 equation (13)], 

(4.10) 0<E^i-Kt^)<c/{l + \nt^), 

we get 

oo oo 

j2\m\\ux)\{XnyEyu-^nt'^n < 5^i/»ii0„(x)|ai 



1 -|- Xnt"^ 

(4.11) < ct--J2\fin)\\Mx)\K7'<oo, 



n=l n=l 

oo 



n=l 



where the last inequality follows from the absolute and uniform convergence of the 
eigenfunction expansion of A™~^/. 

A similar argument using [291 Equation (17)] 



(4.12) 



dt 



x„f- 



1 



and the fact that the eigenfunction expansion of A'"~^/ converges absolutely and 
uniformly allows us to differentiate the series ( 14. Sp term by term with respect to t. 

We next show that A'm G C^{D) for / = l,---,m — 1: this follows from the bounds 
in [211 Theorem 8.33] and the absolute and uniform convergence of the series defining 



(4.13) 



\(t>n\l,a;D < C(l + A„) SUp |(/)„(x) | , 

D 
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where C = C{d, dD) is a finite constant. Here 

k 

\u\k,a;D = sup [D'^ u]a^j:i + sup sup \D"'u\, A; = 0, 1, 2, ■ ■ ■ 

and 

[D^u]^,D= sup p— — 

x,y£D,Xyiy !J\ 

are norms on C^''"(-D). Hence for / = 0, 1, 2, ■ ■ ■ , m — 1 

oo 

|A'n(.,t)|i,,;Z5 < cV|/H|E/3(-A„t'^)(A„)'(l + A„)sup|</.„(x)| 

< C5^|/(n)(A„)'|^^y^sup|</.„(x)| 

n=l 

oo 

< Ct-^Vsup|0„(x)|(A„)V»l 

oo 

+C J]sup|0„(x)|(AJV»l < oo- 



1 D 

n=l 



Witli tliis we establislied tliat u G "Ha- n Cb{D^) fl C^(Z:'). 
Observe next tliat tlie Laplace transform of 

a m-1 ,j/jy^_l 



IS 

l/m "1-1 



1 + V s-^/™f-A V - ^~ 



„l/7n-l 



0, 



since the Laplace transform of A^t^/*") is ^f/J^^ and using (14.31) . Hence, 

the uniqueness of Laplace transforms, we get 



^ T{j/m) 
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Now applying the time derivative and A™ to the series in fl4.8l) term by term gives 



9 < ^ 



n=l 

oo 

n=l 





V ( 

^ r(j7m) 



(x)-Ei 



9_ 
at 



E F777z:T(-^n)' - (-An)'"i?j,(-A.t^) 



^ r(j/m) 



which shows that the PDE in f ll.lOp is satisfied. Thus, we conclude that u defined by 
f l4.8p is a classical solution to fll.lOp . This completes the first part of the proof. 

We next deal with the case m = 2^ Observe that E^/^\t), process inverse to stable 
subordinator of index 1/2^ and have the same density, as they have the same 

one dimensional distribution by Lemma 13.21 Let /) denote the common density 
of and E^/^\t). 

oo 

(4.14) u{t,x) = Y,f{n)Ey^u{-Kt^l^'')cl>r,{x) 



n=l 

"OO 







(4.15) 



(4.16) 



n=l 

oo 

n=l 



n]e 



n e 



-Xnl. 



[X 



p{t,l)dl 



Tn{l)f{x)p{t,l)dl 
E.[/(X(|4(t)|))/(r^(X) > |4(t)|)]. 



Note that Equation f l4.16p follows by a conditioning argument. 

Finally we prove uniqueness. Let Ui, U2 be two solutions of fll.lOp with initial data 
m(0,x) = f{x) and Dirichlet boundary condition u{t,x) = for x G dD. Then 
U = Ui — ^2 is a solution of f ll.lOp with zero initial data and zero boundary value. 
Taking 0,„-transform on both sides of f ll.lOp we get 



d_ 

dt 



U{t, n) = {-\nTU{t, n), [/(O, n) = 0, 



and then U{t,n) = for all t > and all n > 1. This implies that U{t,x) = in 
the sense of functions, since {0„ : n > 1} forms a complete orthonormal basis 
for L'^{D). Hence, U(t,x) = for all t > and almost all x & D. Since [/ is a 
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continuous function on D, we have U{t,x) = for all {t,x) G [0, oo) x D, thereby 
proving uniqueness. □ 

Corollary 4.1. Let f G C^^{D) he a 2k-times continuously differentiable function of 
compact support in D. Ifk>m — l + 3d/4, then the Equation i\1.10\) has a strong 
solution. In particular, if f ^ C^{D), then the classical solution of Equation Iil.l0\) 
is in C^{D). 



Proof. By Example 2.1.8 of [19], |0„(x)| < {XnY'^. Also, from Corollary 6.2.2 of [20 
we have A„ ~ n'^/'^. 
By Lemma [4.11 we get 



(4.17) AV(n) = (-A„)V». 

Using Cauchy-Schwartz inequality and the fact / G C^^{D), we get 



{A'^f{x)fdx 



D 



1/2 r 



{4>n{x)Ydx 



D 



1/2 



D 



1/2 



Cfe, 



where is a constant independent of n. 

This and Equation fliTTI) give \f{n)\ < Cfc(A„)-^ 
Since 

oo 

A'"^V(x) = J2i-Xnr~'f{n)<PM, 

n=l 

to get the absolute and uniform convergence of the series defining A'""^/, we consider 

oo oo 

5^(A„)"^-Vn(x)||/»| < 5^(A„)'^/^+"^-^c,(A„,)-('=) 



n=l 



n=l 



< Cn y^(r^2/d^d/4+m-l-fc = Cfc ^ ^l/2+2(m-l)/d-2A.-/d 



^2/d'jd/4+m-l-fc 
n=l n=l 

which is finite if -2(m - l)/rf - 1/2 + 2k/d > 1, i.e. k > m - 1 + 3d/4 



□ 



5. Other higher order Cauchy Problems 

The (i- dimensional symmetric a-stable process X{t) with a G (0,2] is the process 
with stationary independent increments whose transition density 

p^x, y) = p'^it, x-y), (t, X, y) G (0, oo) x R"' x 

is characterized by 

f e^^V(t,l/)c?l/ = exp(-t|er), t>0,eGR'^. 

The process has right continuous paths, it is rotation and translation invariant. For 
a = 2, this is Brownian motion running twice the speed of standard Brownian motion. 
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Since the Laplace transform method does not apply by the appearance of t ^ in 
the PDE fll.l2p . we use a direct method to prove Theorem II .31 

Proof of Theorem \1.3[ By a simple conditioning argument and using the series rep- 
resentation of the killed semigroup TD{t) in (12. 2p . we can express u{t,x) as 

(5.1) = 2f[±e-^-fM))^^^Y''- 
We use Fubini-Tonelli theorem, the simple inequality 

POO . 1 

(5.2) / e""-- ' ds < 



7r(t2 + s2) - ntXn 

and the fact that the series defining A/ is absolutely and uniformly convergent to 
show that 

°° POO ^ 

x) = 2j2 /»0n(x) / e-""- ds. 
^ Jo 7r(t^ + s') 

Using this, and (15. 2 p we can show that we can apply to the series (15. ip term by 
term. Hence it follows that u, Au, A^n G C{Doo)- 

We next show that each term in the series (15. ip satisfy the PDE (II. 3p . We use the 
fact that {t, s) satisfy 

(5^ + ajj)p'(M) = o. 

dominated convergence theorem, and integration by parts twice to get 

n2 / _ poo \ /-oo o2 

— e-'-V{t,s)ds] = f{n)Ux) e-'-'—p\t,s)ds 

poo o2 

= ^An0„(x)/(n) 

7lt 

POO 



+ A^0„(x)/(n) / e-^-Y{t,s)ds 
1 







, , ,.(2;)/(n) 

Tit 

POO 

+ -AVn(x)/(n) / e-^"y(t,s)rfs. 
Jo 

From this we get that the time derivative can be applied term by term to the series 
(15. ip since the series for A/ and A^m converge absolutely and uniformly. Hence 
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applying the time derivative and the term by term to the series fl5.ip we obtain 



oo 
n=l 



nt 



1 



Ot^ TXt 



(5.3) 



+AVn(x) / 

Jo 



e-^-'p\t,s)ds ) = 0. 



□ 



For a rational a ^ 1 the PDE is more complicated since kernels of symmetric 
a-stable processes satisfy a higher order PDE. 

Theorem 5.1 (Nane [12])- Let a G (0,2) be rational a = l/m, where I and m are 
relatively prime. Let T{s)f{x) = Kx[f{X{s))] be the semigroup of Brownian motion 
X{t) and let A be its generator. Let f be a bounded measurable function in the 
domain of A, with D'^f bounded and Holder continuous for all multi index 7 such 
that I7I = 21. Then u{t,x) = Kx[f{X{\Y{t)\)] is a solution of 

i=l ^ ^ 

- A'^^u{t,x), t>0, xeR'^ 
u{0,x) = fix), xgM°'. 



Let 



]C^2i = \u:Doo^M: -^u, A^u, 6 C(Doo) for j = 1, ■ ■ ■ , 21, 



A^u,eCiD) for J = I,-- - ,2/- 1 



This theorem takes the following form in bounded domains 

Theorem 5.2. Let a G (0,2) be rational a = l/m, where I andm are relatively prime. 
Let D be a bounded domain with dD G C^''^,0 < 7 < 1. Let {X(t)} be a Brownian 
motion, and let {Yit)} be symmetric a-stable process independent of {X(t)}. Let 
f G -D(A) n C^'~^(D) for which the eigenfunction expansion of A^''~^f with respect 
to the complete orthonormal basis {(pn '■ n > 1} converges absolutely and uniformly. 
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Then u{t,x) = Kx[f{X{\Y(t)\)I{T£){X) > \Y{t)\)] is a classical solution of 

i=l ^ ^ 

- A^'u(t,x), t > 0, xe D 
u{0, x) = f{x), X E D. 

u{t,x) = A%(t,x)=0, X edD, t>0, j = !,■■■ ,21-1. 

Proof. Let TD{t) be the killed semigroup of Brownian motion in D. We use the 
representation 

POO 

u{t,x) = E,[f{X{\Y{t)\)I{rD{X)>\Y{t)\)] = 2 p^{t, s)To{s)f{x)ds 

Jo 



(5.5) 



POD / °° \ 

2/ i^J2e~'-Jin)M^)jp'^{t,.s)ds 



and the fact that the transition density s) of the process Y{t) satisfies 

ig^y + = 0, (t, x) e (0, 00) X M, 



from Lemma 3.2 in 

We also use the well-known fact that for a G (0, 2) 

where 

(27r) 

Hence 



This allows us to deduce the fact that A% G C{D^) for j = 0, 1, ■ ■ ■ ,2/ and also 
that A% G C{D) for j = 0, 1, ■ ■ ■ , 2/ - 1. 

To get that ■§^u{t,x) G C{Dao) we use the fact that the series defining A^'u 
converges absolutely and uniformly as well as the series defining A^^~^/ converges 
absolutely and uniformly and the fact that the terms in the series defining u{t,x) in 
(1531) satisfy the PDE CT . 

Hence we can interchange the sum and powers of the Laplacian in the series (15.51) 
term by term to show that PDE in (15.41) is satisfied. □ 
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6. Extensions and Discussion 



A uniformly elliptic operator of divergence form is defined on functions by 

d {aij{x){du/dxi)) 



(6.1) Lu=Y, 



i,j=l 



dx 

with aij{x) = aji{x) and, for some A > 0, 



(6.2) >^J2y'^ E «^^(^)?/^%- < E^'' ^ 

The operator L acts on the Hilbert space L'^{D). We define the initial domain 
C^{D) of the operator as follows. We say that / is in C^{D), if / e C°°{D) and 
f{x) = for all X G dD. This condition incorporates the notion of Dirichlet boundary 
conditions. 

From [20i Corollary 6.1], we have that the associated quadratic form 



is closable on the domain C^{D) and the domain of the closure is independent of the 
particular coefficients (ajj) chosen. In particular, for f,g E C^{D) by integration by 
parts 



g{x)Lf{x)dx = Q{f,g) = / f{x)Lg{x)dx, 

'D J D 

which shows that L is symmetric. 

From now on, we will use the symbol Lo if we particularly want to emphasize the 
choice of Dirichlet boundary conditions, to refer to the self-adjoint operator associated 
with the closure of the quadratic form above by the use of [20| Theorem 4.4.5]. Thus, 
Lj) is the Friedrichs extension of the operator defined initially on C^{D). 

If the coefficients aij{x) are smooth (ajj(x) G C^{D)), then L^u takes the form 

d r^9 d / d r\ / ^ 09 *^ O 

Lnu = ^ a.(x)^+5: E ^ = E «'^(-)^+E ^^(-)^- 

i,j=l * ^ i=l \j=l ■> / ij=l * J j=l * 

If Xt is a solution to 

dXt = aiXt)dWt + 6(X0c?t, Xo = xo, 

where a is a (i x c? matrix, and Wt is a Brownian motion, then Xt is associated with 
the operator with a = cra^, (see Chapters 1 and 5 of Bass [I2])- Define the first 
exit time as td{X) = mi{t > : Xt ^ D}. The semigroup defined by T(t)f{x) = 
Ex[f (Xt) I {t£,{X)) > t)] has generator L^,, which follows by an application of the Ito 
formula. 
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Let D he a bounded domain in M'^. Suppose L is a uniformly elliptic operator of 
divergence form with Dirichlet boundary conditions on D, and that there exists a 
constant A such that for all x G -D, 

d 

(6.3) 5^ |a,,(x)| < A. 

Let T£){t) be the corresponding semigroup. Then Tnit) is an ultracontractive semi- 
group (even an intrinsically ultracontractive semigroup), see Corollary 3.2.8, Theorem 
2.L4, Theorem 4.2.4, and Note 4.6.10 in [12]. Every ultracontractive semigroup has 
a kernel for the killed semigroup on a bounded domain which can be represented as 
a series expansion of the eigenvalues and the eigenf unctions of L^, (cf. [191 Theo- 
rems 2.1.4 and 2.3.6] and [2H Theorems 8.37 and 8.38] ): There exist eigenvalues 
< /ii < /i2 < /is ■ ■ ■ , such that /i„ — ?■ oo, as n — )■ oo, with the corresponding 
complete orthonormal set (in Hq) of eigenfunctions ipn of the operator Ld satisfying 



(6.4) Lotpnix) = -lXntpn{x), X E D : tpnldD = 0. 

In this case, 

oo 

PD{t,x,y) = ^e-'^"Vn(a;)^„,(|/) 

n=l 

is the heat kernel of the killed semigroup T^. The series converges absolutely and 
uniformly on [to, C)o) x D x D for all to > 0. 

Suppose D satisfies a uniform exterior cone condition. Let {Xt} be a Markov 
process in with generator Lo, and / be continuous on D. Then the semigroup 

Tnit)f{x) = E,[f{Xt)I{t < ToiX))] 



(6.5) Jd 



PD{t,x,y)f{y)dy 
5^e-^"Vn(x)/(n) 



n=l 

solves the Dirichlet initial-boundary value problem in D: 

= LDu{t,x), xeD,t>0, 

en 

u{t, x) = 0, X e dD, 
u{0,x) = f{x), xeD. 

Remark 6.1. The eigenfunctions belong to L°°{D) fl C'^{D) for some a > 0, by 
[211 Theorems 8.15 and 8.24]. If D satisfies a uniform exterior cone condition all the 
eigenfunctions belong to C^iD) by [21 Theorem 8.29]. Ifaij e C^^D) anddD e C^'", 
then all the eigenfunctions belong to C^'"(Z)) by [211 Corollary 8.36]. If Qij G C°°{D) 
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then each eigenfunction of L is in C°°{D) by [211 Corollary 8.11]. If aij G C°°{D) 
and dD G C°°, then each eigenfunction of L is in C°°{D) by pH Theorem 8.13]. 

Remark 6.2. Theorems \l.l\ li.jjl and \l.'^ are valid if we replace the outer Brownian 
motion with a Diffusion process. This can be verified by considering fl4.13p . fl6.4p . 
(16.51) . Remark \6.1\ and Theorem 3.1 in [37] . 

Remark 6.3. It might be an interesting project to consider the PDEs treated in this 
paper with the Neumann boundary conditions. Probably the solutions will be obtained 
by running reflected diffusions subordinated by k-iterated Brownian motions. We will 
treat this problem elsewhere. 

Acknowledgments. I would like to thank Mark Meerschaert for encouragement 
on working on the problems in this paper. I also would like to thank anonymous 
referee for his or her helpful comments on the comments on the historical priority on 
the equivalence of fractional Cauchy problems and higher order PDEs which improved 
the accuracy of the results leading to the present paper. 

References 

[1] Agrawal, O.P. (2002). Solution for a fractfonal diffusfon-wavc equation defined in a bounded 

domain. Fractional order calculus and its applications. Nonlinear Dynam. 29 145-155. 
[2] AUouba, H., A linearized Kuramoto-Sivashinsky PDE via an imaginary- Brownian-time- 

Brownian- angle process, C.R. Acad. Sci. Paris Ser. 1336 (2003), 309-314. 
[3] AUouba, H., Brownian-time processes: The pde connection and the corresponding Feynman-Kac 

formula, Trans. Amer. Math. Soc. 354 (2002), no. 11 4627 - 4637. 
[4] AUouba, H. and Zheng, W. (2001). Brownian-time processes: The PDE connection and the 

half-derivative generator, Ann. Prob. 29 1780-1795. 
[5] Applebaum, D. (2004). Levy Processes and Stochastic Calculus. Cambridge studies in advanced 

mathematics. 

[6] Arcones, M. A. The large deviation principle for stochastic processes. I. Teor. Veroyatnost. i 
Primenen. 47 (2002), no. 4, 727-746; translation in Theory Probab. Appl. 47 (2003), no. 4, 
567-583. 

[7] Arendt, W., Batty, C, Hieber, M. and Neubrander, F. (2001). Vector-valued Laplace transforms 

and Cauchy problems. Monographs in Mathematics, Birkhauser-Verlag, Berlin. 
[8] Aurzada, F. and Lifshits, M. (2009). On the Small Deviation Problem for Some Iterated Pro- 
cesses. Electronic. J. Probab. Vol. 14, Paper no. 68, pages 19922010. 
[9] Baeumer, B. and Meerschaert, M.M. (2001). Stochastic solutions for fractional Cauchy prob- 
lems. Fractional Calculus Appl. Anal. 4 481-500. 
[10] Baeumer, B. Meerschaert, M.M. and Nane, E. Brownian subordinators and fractional Cauchy 

problems. Trans. Amer. Math. Soc. 361 (2009), 3915-3930. 
[11] Bass, R.F. (1995). Probabilistic Techniques in Analysis. Springer- Verlag, New York. 
[12] Bass, R. F. (1998). Diffusions and Elliptic Operators. Springer- Verlag, New York. 
[13] Bcckcr-Kern, P., Meerschaert, M.M. and Scheffler, H.P. (2004). Limit theorem for continuous 

time random walks with two time scales. J. Applied Probab. 41 No. 2, 455-466. 
[14] Benachour, S.; Roynette, B.; Vallois, P. Explicit solutions of some fourth order partial differ- 
ential equations via iterated Brownian motion. Seminar on Stochastic Analysis, Random Fields 
and Applications (Ascona, 1996), 39-61, Progr. Probab., 45, Birkhuser, Basel, 1999. 

24 



[15] Bochncr, S., Harmonic Analysis and the Theory of Probability, Univ. California press, Berkeley. 
CA (1955). 

[16] Burdzy, K. (1993). Some path properties of iterated Brownian motion. In Seminar on Stochastic 
Processes (E. Qinlar, K.L. Chung and M.J. Sharpe, eds.) 67-87. Birkhauser, Boston. 

[17] Burdzy, K. (1994). Variation of iterated Brownian motion. In Workshops and Conference on 
Measure-valued Processes, Stochastic Partial Differential Equations and Interacting Particle 
Systems (D.A. Dawson, ed.) 35-53. Amer. Math. Soc. Providence, RI. 

[18] Caputo, M. (1967). Linear models of dissipation whose Q is almost frequency independent. Part 
II. Geophys. J. R. Astr. Soc. 13 529-539. 

[19] Davies, E.B. (1989). Heat Kernels and Spectral Theory. Cambridge Univ. Press, Cambridge. 

[20] Davies, E.B. (1996). Spectral Theory and Differential Operators. Cambridge Univ. Press, Cam- 
bridge. 

[21] DcBlassie, R.D. (2004). Iterated Brownian motion in an open set. Ann. Appl. Probab. 14 1529- 
1558. 

[22] DiBcnedetto, E. (1995). Partial Differential Equations. Birkhauser, Boston. 

[23] Funaki, T. A probabilistic construction of the solution of some higher order parabolic differential 

equations, Proc. Japan Acad. Ser. A. Math. Sci. (5) 55 (1979), 176 - 179. 
[24] Gilbarg, D. and N.S. Trudinger, N.S. (2001) Elliptic Partial Differential Equations of Second 

Order. Reprint of the 1998 ed.- New York Springer. 
[25] Gorenflo, R. and Mainardi, F. (2003). Fractional diffusion processes: Probability distribution 

and continuous time random walk. Lecture Notes in Physics 621 148-166. 
[26] Hochberg, K. J. and Orsingher, E. (1996). Composition of stochastic processes governed by 

higher-order parabolic and hyperbolic equations. J. Theoret. Probab. 9 511532. 
[27] Jacob, N. (1996) Pseudo- Differential Operators and Markov Processes. Berlin : Akad. Verl. 
[28] Khoshnevisan, D. and Lewis, T.M., Chung's law of the iterated logarithm for iterated Brownian 

motion, Ann. Inst. H. Poincare Probab. Statist. 32 (1996), no. 3, 349-359. 
[29] Krageloh, A.M. (2003). Two families of functions related to the fractional powers of generators 

of strongly continuous contraction semigroups, J. Math. Anal. Appl. 283 459-467. 
[30] Le Gall, J.-F. (1993). Solutions positives de An = dans le disque unite. C.R. Acad. Sci. Paris 

Ser. I 317 873-878. 

[31] Le Gall, J.-F. (1994). A path-valued Markov process and its connections with partial differential 

equations. In Eurepean Congress of Mathematics 2 185-212. Birkhauser, Boston. 
[32] Le Gall, J.-F. (1995). The Brownian snake and solutions of Aw = ii^ in a domain. Probab. 

Theory Related Fields 102 393-432. 
[33] Mainardi, F. and Gorenflo, R. (2000). On Mittag-Lcffer-type functions in fractional evolution 

processes, J. Comput. Appl. Math. 118 283-299. 
[34] Mecrschaert, M.M. and Scheffler, H.P. (2001) Limit Distributions for Sums of Independent 

Random Vectors: Heavy Tails in Theory and Practice. Wiley Interscience, New York. 
[35] Mecrschaert, M.M. and Scheffler, H.P. (2004). Limit theorems for continuous time random walks 

with infinite mean waiting times. J. Applied Probab. 41 No. 3, 623-638. 
[36] Mecrschaert, M.M., Benson, D.A., Scheffler, H.P. and Baeumer, B. (2002). Stochastic solution 

of space-time fractional diffusion equations. Phys. Rev. E 65, 1103-1106. 
[37] Mecrschaert, M.M., Nane, E. and Vellaisamy, P. (2009). Fractional Cauchy problems in bounded 

domains, Ann. Prob. Vol. 37, No. 3, 979-1007 
[38] Mecrschaert, M.M., Nane, E. and Xiao, Y., Correlated continuous time random walks, Stat. & 

Probab. Lett. Volume 79, Issue 9, 1 May 2009, Pages 1194-1202 
[39] Nane, E., Iterated Brownian motion in parabola- shaped domains, Potential Analysis, 24 (2006), 

105-123. 



25 



Nanc, E., Iterated Brownian motion in bounded domains in R", Stochastic Processes and Their 
Apphcations, 116 (2006), 905-916. 

Nane, E., Lifetime asymptotics of iterated Brownian motion in R": Esaini:P&S, March 2007, 
Vol. 11, pp. 147-160. 

Nanc, E., Higher order PDE's and iterated processes, Trans. Amer. Math. See. 360 (2008), 
2681-2692. 

Nigmatuhin, R.R. (1986). The reahzation of the gencrahzed transfer in a medium with fractal 
geometry. Phys. Status Solidi B 133 425-430. 

Nourdin, I. and Peccati, G., Weighted power variations of iterated Brownian motion. Electronic 
J. Probab. Vol. 13 (2008), 1229-1256. 

Orsingher, E. and Beghin, L., (2004) Time-fractional telegraph equations and telegraph pro- 
cesses with Brownian time. Proh. Theory Rel. Fields 128, 141-160. 

Orsingher, E. and Beghin, L., (2008) Fractional diffusion equations and processes with randomly 
varying time, Ann. Probab. Volume 37, Number 1 (2009), 206-249. 
Podlubny, I. (1999). Fractional Differential Equations, Academic Press, San Diego. 
Royden H.L. (1968). Real Analysis. 2nd Edition, MacMillan, New York. 

Samko, S., Kilbas, A. and Marichev, O. (1993). Fractional Integrals and derivatives: Theory 
and Applications. Gordon and Breach, London. 

Sato, K.I. (1999). Levy Processes and Infinitely Divisible Distributions. Cambridge University 
Press. 

Scalas, E. (2004). Five years of Continuous- Time Random Walks in Econophysics. Proceedings 
of WEHIA 2004, A. Namatame (ed.), Kyoto. 

Schilling, R.L., Growth and Holder conditions for sample paths of Feller proceses. Probability 
Theory and Related Fields 112, 565-611 (1998) 

Schneider, W. R. and Wyss, W. (1989). Fractional diffusion and wave equations. J. Math. Phys. 
30, 134-144. 

Song, R. and Vondracek, Z. (2003). Potential theory of subordinate killed Brownian motion in 
a domain. Probab. Theory Relat. Fields 125 578-592. 

Zaslavsky, G. (1994). Fractional kinetic equation for Hamiltonian chaos. Chaotic advcction, 
tracer dynamics and turbulent dispersion. Phys. D 76 110-122. 

Erkan Nane, Department of Mathematics and Statistics, 221 Parker Hall, Auburn 
University, Auburn, AL 36849 
E-mail address: nane@auburn.edu 



26 



